Generation of ring dark solitons by phase engineering and their oscillations in spin-1 

Bose-Einstein condensates 
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The ring dark solitons in spin-1 ^'^Na and *^Rb Bose-Einstein condensates are studied numeri- 
cally in the framework of the time-dependent Gross-Pitaevskii equations. By simulating the phase 
engineering technique in real experiments, we explore the roles of the parameters characterizing the 
far-off resonant laser pulse which can be used to generate the ring dark solitons. The variations 
of these parameters have dramatic effect on the lifetime and the decay profiles of the ring dark 
solitons. If only one ring dark soliton is generated in one component of the condensate, ring dark 
solitons in other components are inclined to be induced, resulting in a coexistence state composed 
of interdependent ring dark solitons coming from different components of the condensate. Ring 
dark solitons in this coexistence state exhibit dynamical oscillations for hundreds of milliseconds. 
By studying the lifetime and decaying profiles of ring dark solitons, we explore the similarities and 
differences of ^^Na and *^Rb condensates. Besides, taking into account the fact that the center of 
the ring may not be coincide with that of the trap, we study the dynamics and decaying profiles of 
the off-centered ring dark solitons in the presence of symmetry breaking effect. 

PACS numbers: 03.75.Lm, 03.75.Mn, 05.45.Yv 



I. INTRODUCTION 

Since the experimental realization of the Bose-Einstein 
condensate (BEC), intensive interest has been growing 
and the condensate has been a popularly investigated 
platform for various effects of quantum many-body inter- 
action. This opens a new field of exploring various types 
of excitations such as bright solitons, dark solitons and 
vortices [1141. The concept of ring dark soliton was intro- 
duced theoretically and realized experimentally in non- 
linear optics 0-01 • Theocharis et al. investigated the ring 
dark solitons in scalar BEC with disk-shaped trap 04111 ■ 
Recently, dark-bright stripe solitons and ring dark soli- 
tons in two-component BEC are studied and the stable 
ring dark-bright solitons are found [l^, [I!]- The deep 
ring dark soliton in scalar BEC has been found always 
unstable towards the formation of vortices. However, it 
is shown in [l^] that the presence of the "bright" compo- 
nent has a stabilizing effect on the symbiotic state. 

As far as we know, ring dark solitons in spin-1 BEC 
have not been investigated theoretically or experimen- 
tally yet. Spinor BEC liberates the spin degrees of free- 
dom of the atoms confined in optical traps. Several 
phases are possible below the transition temperature Tc, 
and this kind of diversity leads to a fascinating area of 
research for quantum gases. The magnetic moment ac- 
companying the spin gives rise to magnetism, resulting in 
kinds of excitations like spin domains, spin textures and 
fractional vortices ■ In addition, matter- wave dark 

solitons can be created experimentally by means of vari- 
ous methods, including the phase engineering and density 
engineering techniques |16l - [2(| . Especially, Wright et al. 



[2ll | have achieved spatially varying control of the ampli- 
tude and phase of the spinor order parameter. To shade 
light on the manipulation of ring dark solitons in spinor 
BEC, we investigate the generations and oscillations of 
the ring dark solitons in spin-1 BEC by simulating the 
phase engineering technique. 

In the present paper, we explore the dynamical genera- 
tion and evolution of the ring dark solitons in = 1 ^'^Na 
and ^^Rb condensates. By using models to simulate the 
phase engineering technique in the experiments, we study 
the effect of two key parameters of the far-off resonant 
laser pulse on the dynamics of the ring dark solitons, the 
width of the laser-induced Stark-shift potential and the 
duration of the laser beam. We study the generation and 
dynamics of the coexistence state of ring dark solitons 
that is induced by producing only one ring dark soliton 
in one component of the condensate. Numerical results 
are compared between ^'^Na and ®^Rb condensates. 

The rest of the paper is organized as follows. In Sec. 
im we introduce the mean field theoretical model of spin- 
1 BEC, and a phase distribution model to imitate the 
phase configuration of the order parameter after the il- 
lumination of the laser beam. In Sec. IIIIl we obtain 
the ground states of the spin-1 ^^^Na and ^^Rb conden- 
sates numerically, then study the generation process of 
the ring dark solitons and investigate the roles of param- 
eters characterizing the far-off resonant laser pulse. Sec. 
IIVI focuses on the emergence, oscillation and decay of the 
coexistence state of ring dark solitons that is induced by 
the ring dark soliton in one component of the condensate. 
Finally, the conclusions are summarized in Sec. [V] 



FIG. 1: (Color online) A schematic of imprinting phase steps 
onto the condensate. The far-off resonant laser pulse (light 
red) goes through an absorption plate (gray), which is used 
to modulate the laser. The modulated laser pulse illuminates 
the cold atomic condensate (blue) and produces an expected 
phase configuration. Here, Ro is the radius of the ring dark 
soliton that is to be created. 



FIG. 2: (Color online) (a) The ID density distribution 
lipiix^y — 0)p showing the emergence of a ring dark soli- 
ton in the tjji component of the ^''Na condensate with a total 
number of atoms Nnu = fO^, Ax = 0.54 fim, A(f> = n and 
-Ro = 15.77 /xm. (b) The ring dark soliton configuration in 
the ipi component of the ^^Na condensate at t = 0.64 ms. For 
the ^'^Rb condensate, the process of generating and configu- 
ration of the ring dark soliton show no significant difi'erence 
comparing to the ^^Na condensate. 



II. THE THEORETICAL MODEL 

Most BECs of atomic gases have internal degrees of 
freedom originating from the spin. For the *^Rb and 
^■^Na atoms, there are spin combinations with S = 1/2 
(electronic spin) and 1 = 3/2 (nuclear spin), so the to- 
tal spin F can be chosen to be = 1 or = 2. Both 
spin-1 and spin-2 atomic condensates are investigated ex- 
perimentally and theoretically [3, [22| - [3l| . For F = 1, 
there arc three sub-level states with magnetic quantum 
number m = 0,±1. Within the mean field description, 
the spin-1 cold atom condensate can be described by a 
vectorial order parameter ^ = {■ipi,ipo,ip-i)'^ (the su- 
perscript T stands for the transpose), where the com- 
ponent ipm (jn = 0, ±1) denotes the macroscopic wave 
function of the atoms condensed in the hyperfine state 
\F = 1,to). The dynamics of the vectorial order param- 
eter follows the multicomponent Gross-Pitaevskii equa- 
tions (GPEs)!!? 
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where V2 = d^/dx^ + d^/dy^ + d'^/dz^, co = 4:TTff{aQ + 
2a2)/3M and C2 = 47r/i^(a2 — ao)/3Af (M is the atomic 
mass), ao and 02 are the s-wave scattering lengths cor- 
responding to total spin of the two colliding bosons 
and 2, respectively. The ground state of the conden- 
sate is ferromagnetic for C2 < and antiferromagnetic 
for C2 > 0. V = M{ujIx^ +ujly^ +ujIz^)I2 [u^.y,^ is the 
confinement frequency in the corresponding direction) is 
the external potential. The density of particles is defined 



n = X^rjj ~ 0, ±1) and the total particles 

number N ~ J ndr. 

To investigate the ring dark solitons in a pancake con- 
densate, we choose lOz (jJx = ^y, and the 3D GPEs 
([1]) can be reduced to 2D GPEs. After the dimensionless 
process (we denote the order parameter with the same 
■0m), we have 
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where a,. 



N^/WiJJ(2Tr¥)co, as 



N y/MhJ^/(2TTW)c2 and il = uix/i^z- The particle num- 
ber conservation condition becomes / X)m IV'mP'^'^ = 1 
(to = 0, ±1). The time and length arc measured in units 
of and y/h/ (MuJz), respectively. Based on the 2D 
GPEs ([2]), we are going to shade light on the dynamical 
properties of ring dark solitons in spin-1 BEG. 

Experimentally, the matter wave dark solitons can be 
created by means of various methods, such as phase en- 
gineering, density en gine ering and quantum-state engi- 
neering 0, [la - llQl I21I |32| . Here we study the dynamics 
of ring dark solitons by simulating the phase engineering 
technique |1, El [13, El . 

We suppose that ujj_ ~ lo^ = ujy = (27r) x 12 Hz and 
u)z — (27r) X 372 Hz, respectively. As in [33, we choose 
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FIG. 3: (Color online) The ID density |V'i(a:,y = 0)^ (solid 
line) and phase distribution (^(x, y = 0) (dotted line) sX t = 
0.21 ms in the *^Rb condensate with Ni^y = 2.5 x lO"*, Aa; = 
0.54 /im, Ro = 15.77 nm and (a) = 7r/2, (b) = tt, (c) 
A(j!> = 37r/2, respectively. 



ao = 46aB, a2 ~ 52aB (as is the Bohr radius) for ^'^Na 
condensate and oq = llOas, 02 = lOTas for ^""Rb con- 
densate. And we numerically trace the dynamics of the 
ring dark solitons in spin-1 ^'^Na and ^''Rb condensates. 

Generally, the phase engineering technique consists of 
passing a far-off resonant laser pulse through an absorp- 
tion plate which is used to modulate the spacial distribu- 
tion of the intensity of laser beam, and creating the cor- 
responding conservative Stark-shift potential that leads 
to a space-dependent phase shift in the condensate order 
parameter. Fig. [Ijshows the schematic for the generation 
of ring dark solitons by phase engineering technique. The 
laser beam that is applied onto the condensate supplies 
an effective potential, which would produce phase dis- 
tributions in the condensate order parameter. The laser 
beam is pulsed on for a time St, which is much smaller 
than the correlation time of the condensate. The phase 
factor imprinted on the condensate can be adjusted by 
varying St. We use the following model to simulate the 
phase distribution after illustration of the laser beam, 

'\>{^. 2/) = ^ X [1 + tanh( ^ (3) 

where Ac/i, i?o are the total phase difference and the ra- 
dius of ring dark soliton, respectively. Aa; reveals the 
width of potential edge (induced by the laser beam). 



III. STRUCTURES OF RING DARK SOLITONS 

The ground state is antiferromagnetic for the ^^Na con- 
densate and ferromagnetic for the ^''Rb condensate. Be- 
cause all spinors related to each other by gauge transfor- 
mation e'^ and spin rotations Via, /3, 7) (a, /3, 7 are Euler 
angles) are degenerate 1331 , the ground vectorial order pa- 
rameter ($1, <i>Oj in our numerical results is only 
one of them. The ground state of the ^^Na condensate 
has not the component with magnetic quantum number 
TO = 0. In certain cases we use ground state of the ^'^Na 
condensate that is produced by rotating the numerically 
obtained ground state in the spin-1 space. We choose the 
rotating operation as ?7(0, 7r/2, 7r/4) to produce number 
balanced states in $1 and $_i components. 
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FIG. 4: (Color online) The ID phase distribution 0(a;, y = 0) 
at f = 0.43 ms in the t/ji component of the *^Rb condensate 
with parameters Tio = 15.77 pm, A^ = n, Ax = 1.09 fim 
(dashed line) and Aa: = 0.25 /xm (solid line), (a) and (b) show 
the corresponding different density decay profiles of the ring 
dark solitons ai t = 12.84 ms and t — 57.76 ms, respectively. 

After the illustration of the far-off resonant laser pulse, 
the matter wave function can be written as 

* - (Vi>o,V'-if = {^i,<fo,^-ife'^(--^y\ (4) 

where (j){x, y) is the phase distribution after the illumi- 
nation of the laser beam as modeled in equation ([3|). 

As a far-off resonant laser pulse passes through the 
condensate, all the three components of the condensate 
feel the same conservative Stark-shift potential, and ring 
dark solitons in different components come up with the 
same configuration. So we just show one component to 
exhibit the early stages of the dynamics of the conden- 
sate. Using initial conditions 4* in equation ^ with the 
numerical ground state before the rotation in the spin 
space, we show the generation process of the ring dark 
solitons in the i/'i component of ^'^Na condensate in Fig. 
[2] with a total number of atoms A^jva = 10^, A0 = tt. 
Ax = 0.54 and i?o = 15.77 /xto. As the total phase 
difference is imprinted onto the condensate, the density 
in each component begins to redistribute, and one obvi- 
ous ring dark soliton appears, as shown in Fig. [5] (a). 
At about t = 0.02 ms, the ring dark soliton begins to 
appear, and a bright one is forming up simultaneously 
inside the dark ring. Actually, the relative position of 
the dark and bright ring is related to the phase gradient 
properties, as found in 0] . 

There are two key parameters in the phase engineering 
technique, Ax which reveals the width of the potential 
edge and the phase difference A;/). As the numerical re- 
sults in Fig. 13] show, the phase difference at < = 0.21 ms 
is smaller than the corresponding initial value. In fact, 
the phase gradient increases the total energy of the sys- 
tem, and the depleted phase differences are compensated 
in the form of the fringes of the density distribution as we 
can see from the Fig. [3] (b) and (c). When A(/) = 37r/2, 
a second much shallower ring dark soliton appears, as 
shown in Fig. [3](c). For (a) and (b) in Fig. [3l the parti- 
cle density in the trough chooses a finite value (meaning 
not a black one), whereas the generated ring dark soli- 
ton takes a nearly black form with respect to the density 
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our system, ^ « 0.55 /xm, so it is reasonable to assume 
the unstability of the soliton structures with Ax < 0.55 
ms. In Fig. 21 the notable oscillation of the phase distri- 
butions around the slope enhances the decay of the ring 
structures. 
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FIG. 6: The ID density (top panels) and phase (bottom pan- 
els) distribution ai t — 0.43 ms for the three components of 
the **''Rb condensate, with N^b = 2.5 x 10*, Ro = 15.77 ^Im, 
A4> = 37r/2 and Aa; = 1.09 ^m. 



FIG. 5: (Color online) The density distribution IV'iP of the 
decaying ring dark solitons in the ^^Na and *^Rb condensates 
with parameters N]\!a = A''^;, = 2.5 x 10'*, Ro ~ 8.16 fim, 
A(/!> = 1.057r and Aa; = 0.56 fj.m. Other components of the 
condensate exhibit the same dynamic configuration. The ini- 
tial condition is modeled by equation Q with the numerical 
ground state before rotation in the spin space. 



characteristics in Fig. [3] (c). When we further increase 
the total phase difference, such as Acfi = 2tt, Stt, more 
ring dark solitons will come up and show collective oscil- 
lations. 

Next, we aim to find out the effects of the Stark-shift 
potential edge width on the dynamics of the ring dark 
solions through varying Ax. As Ax varies from 2.17 /im 
to 0.25 /im, the generated ring soliton gets darker and 
darker. But when Ax < 0.54 fim, the lifetime of ring 
dark soliton decreases fiercely and the ring dark soliton 
shows different decay configurations. Using initial con- 
dition modeled in equation (U) for the ^^Rb condensate. 
Fig. Ulshows the phase of the ■01 component in x direction 
with the dashed (solid) line corresponding to Ax = 1.09 
/im (0.25 fj,m). We can see that the phase jump is nearly 
the same, except for small amplitude oscillations for the 
Ax = 0.25 fim case. These two cases correspond to differ- 
ent decay properties of ring dark solions. For Ax — 0.25 
fim, the lifetime of the ring dark solions is shorter, and 
eight vortex and anti-vortex pairs are produced. But for 
Ax — 1.09 fim, the ring dark soliton can last as long as 
57.76 ms, and finally decays into four vortex and anti- 
vortex pairs. This can be understood by resorting to the 
healing length ^ = 1/ \/8nan [lH . The size of topological 
excitations such as dark solitons and vortices is fixed by 
the healing length [HI. And for the *^Rb condensate in 



TABLE I: lifetime of the ring dark solitons in different BEG 
systems. For the scalar BEG with Feshbach resonance man- 
agement (FRM), the total phase difference is A(j> = 0.7077r, 
which should be a key factor for long-lived ring dark soliton. 
For the spinor '^^Na or *^Rb BEG case, ring dark solitons with 
total phase difference Acp « tt are generated in each compo- 
nent of the condensate. 

Systems of BEG Lifetime (units of 1/i^x) 

Scalar BEG without FRM 1.14 
Scalar BEG with FRM 5.09 
Spinor ^^Na BEG without FRM 3.48 
Spinor **^Rb BEG without FRM 4.30 



To explore the differences of the dynamics of the ring 
dark solitons between ^^Na and *^Rb condensates, we 
plot the density distributions of the decaying ring dark 
solitons. As shown in Fig. [5l the density distribution of 
the ipi component of ^^Na and *^Rb condensates exhibits 
the different decaying dynamics. At t = 34.22 ms, nearly 
black ring dark solitons have formed in both condensates. 
However, ring dark solitons in the ^■^Na, which break into 
x-like vortex pairs, live a shorter lifetime comparing to 
that in the ^^Rb condensate which break into -I— like vor- 
tex pairs. The dynamics of the vortex pairs is the same 
as found by in scalar BEC. The configuration of the 
vortex and anti- vortex pairs oscillates between the x- and 
-I— like configurations, while the ring periodically shrinks 
and expands. 

If all the components of the condensate bear the same 
parameterized ring dark solitons, the total density dis- 
tribution of the condensate exhibits the same profile as 
that in each component. When detected in the exper- 
iment with time of flight technique, the profiles of ring 
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FIG. 7: (Color online) The ID density distribution at (a) 
t = 67.60 ms (b) t = 121.50 ms and (c) t = 159.15 ms in 
ipi (solid), and ip-i (dash-dotted) components of the ^"^Na 
condensate showing the "darkness-transferring" effect with 
NNa = 2.5 X 10'*, Ro = 8.16 /im, Afji = 1.057r and Ax = 0.56 
fim. The 2D density distributions according to each time 
point are attached for xpi (left) and ijj-i (right). The density 
of ipo component (which is of the order of 10~*) is omitted. 

dark solitons do not show differences from that in scalar 
BEC. In Tabic U we summarize the lifetimes of ring dark 
solitons in different systems of BEC according to refer- 
ences 0, and our numerical calculations. We note 
that with certain parameters (as shown in Fig. [SJ in ^^Rb 
condensate, the lifetime of the nearly black ring soliton 
in spinor BEC without Feshbach resonance management 
(FRM) is remarkably prolonged referring to the lifetime 
of ring dark soliton in scalar BEC. And this lifetime scale 
is comparable with the lifetime of a shallower ring dark 
soliton (the total phase difference is A(j> = O.TOTtt) in 
[lOj . It would be useful to explore ways to prolong the 
lifetime of the ring dark solitons by varying the parame- 
ters of the external potentials or the spin-dependent and 
independent interactions. 

IV. INDUCED AND OFF-CENTERED RING 
DARK SOLITONS 

A. Common characteristics of the coexistence state 
of ring dark sohtons in ^''Na and *^Rb condensates 

In real experiments, the far-off resonant laser pulse is 
not ideal, so we consider symmetry breaking effects in 
our simulation by allowing the trap center to randomly 
jump within a region [—S, S] x [—S, S] {5 = 0.002/i, where 
h is the grid size). 

We suppose that the condensate takes phase patterns 

as 

* = (V'i,V'o,V'-if = ($ie'^(^'^),<i>o,1>-i)^. (5) 

On this occasion, the ring dark soliton is formed in the 
first component, and at the same time, ring dark solitons 



FIG. 8: (Color online) The time evolution of the radius of the 
ring dark solitons in ^'^Na condensate with iVjva ~ 2.5 x 10*, 
i?o = 8.16 nm and Ax = 0.56 nm. Line a: radius of ring 
dark solitons produced by generating same parameterized ring 
dark solitons in each component of the condensate with an 
initial total phase difference At/i — O.bir; line b: radius of ring 
dark solitons produced by generating one ring dark soliton in 
one component of the condensate with an initial total phase 
difference A(f> — 1.057r. 

in other components are induced. And the mutual fill- 
ing of different components supplies effective potentials 
to support the ring dark soliton structure in each compo- 
nent, resulting in a coexistence structure which can live 
hundreds of milliseconds. 

In Fig. ini we plot the ID density and phase distri- 
butions of the three components at t = 0.43 ms for the 
^""Rb condensate. We find that the ring dark soliton in 
the ipi component has been formed while another two 
much shallower dark solitons in the i/'o and tp-i compo- 
nents germinate as slight phase jumps are in their way 
(see the bottom panel). In fact, the atoms in the ■00 and 
Tp-i components fill the density gap in the ^pl component, 
resulting in a complex interdependent structure. After 
the full generation of ring dark solitons in all the three 
components, the coexistence structure oscillates compli- 
catedly. 

In order to exhibit some aspects of the complex oscil- 
lating properties, we depict the darkness of the solitons in 
different components (the velocity of the dark soliton is 
connected to the darkness) [s^ . In our numerical calcula- 
tion, both ^■^Na and ^^Rb condensates exhibit a darkness- 
transferring property between different components. Fig. 
[7] shows the dynamical process of the transferring for the 
^■^Na condensate using an initial condition modeled by 
equation ^ with ground state that directly obtained by 
imaginary time evolution. At t = 67.60 ms the nearly 
black ring soliton appears in the component. And 
by t = 121.50 ms, the darkness of the ring soliton in the 
ijj-i component has totally been transferred onto the 
component (as shown in Fig. [7](b)). After 37.65 ms the 
ring dark soliton restores its darkness in the 4'i compo- 
nent as shown in Fig. [7] (c) . For the phase distribution 
model in equation ([5]) and parameters in Fig. [7J the co- 
existence structure in both ^^^Na and ^^Rb condensates 
can live more than 800 ms. 

Even though complicated, the oscillations exhibit peri- 
odicity roughly for an initial condition with the absence 
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of the m = component. And the oscihation frequency 
is approximately 0.58W_\/2, which is different from that 
in the scalar BEC 0, 113] ■ The periodical signature can 
be seen from Fig. [51 showing the evolutionary radius of 
the ring dark soliton. In Fig. 8, line a corresponds to the 
state of generating same ring dark solitons in all the three 
components of the condensate, while line b the state of 
generating only one obvious ring dark soliton in one com- 
ponent of the condensate. For line a, the oscillation fre- 
quency is (jJo ~ 2.UQ/\^ (roughly consistent with that 
in for the scalar BEC) while it is Wo « 0.5817/^2 for 
line 5, indicating a prolonged oscillation period. In con- 
trast with the ID condensate, ring dark soliton oscillates 
between the maximum and minimum radius with nearly 
a two times bigger frequency in the 2D condensate. For 
the coexistence state composed of interdependent ring 
dark solitons coming from different components of the 
condensate, the oscillation exhibit periodic character in 
a rough way because of the interaction and radiation of 
the ring dark solitons. During the rough oscillation pro- 
cess, a prolonged oscillation period is found. In fact, the 
oscillation of the ring dark solitons is hindered by the 
embedded particles of another component, resulting in a 
low oscillation frequency of the ring dark solitons. We 
should mention that if we use initial condition with the 
appearance of all the three components of the conden- 
sate, the oscillation is more complicated and so is the 
periodicity. 

If the coexistence state is induced by an obviously off- 
centered ring dark soliton, it suffers an earlier collapse 
for both ^■^Na and ^""Rb condensates with a time scale of 
about 150 ms. In the following, taking into account the 
fact that the center of the ring may not be coincide with 
that of the trap, we suppose an off-centered ring dark 
soliton with an deviation Sp. The oscillation of the ring 
dark solitons are sensitive to such a deviation. When the 
deviation increases to the value Sp =0.20 ms, the lifetime 
of the ring dark solitons decreases evidently. In Fig. ^ 
we show the decaying process of the ring dark solitons 
with 5p = 0.60 fim for ^'^Na. The parameters are the 
same as that in Fig. [TT] (a)-(c). Without deviation Sp, 
the coexistence state of ring dark solitons can live as long 
as 320.88 ms and breaks into four pieces in the ipo compo- 
nent at t = 378.64 ms without the formulation of vortex 
and anti-vortex pairs as shown in Fig. [Il](a)-(c). With 
a finite Sp the coexistence state suffers breakdown at a 
much earlier time. At t = 146.32 ms, the darkness of the 
rings loose their balance, and the darker part of the ring 
is more inclined to decay into vortices. At t = 170.28 ms, 
the ring structure in the ipQ component distorts slightly 
into a form of trapezoid and four vortex and anti- vortex 
pairs come up , as shown by the corresponding top-right 
phase profile panel. And at i = 285.50 ms, the ring dark 
solitons break down totally. The numerical results show 
that the lifetime and dynamics of the ring dark solitons 
are sensitive to such an deviation, indicating that it is 
very essential in the real experiment to control the rela- 
tive position of the optical potential and the far-off reso- 
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FIG. 9: (Color online) The decaying process shown by the 
density and phase distributions of the '^'^Na condensate with 
an initial off-centered ring dark soliton in the component. 
The deviation in y direction Sp is 0.6 and Rq = 8.16 ^im, 
A(j> = Stt/2, Ax = 0.56 fim. At t = 146.32 ms, imbalance in 
the ring dark solitons appears, and 24.96 ms later, the ring 
dark soliton distorts in the -00 component. The ring dark 
solitons are totally broken with vortex and anti-vortex pairs 
in the -i/jo component a.t t — 188.25 ms. 



nant laser beam which is used to generate the ring dark 
solitons. 



B. Different aspects of the coexistence state of ring 
dark solitons in ^'^Na and '^^Rb condensates 

To explore the differences between ^'^Na and ^''Rb con- 
densates in a short time scale, we suppose the initial vec- 
tor parameter as 



(6) 



with the parameters N^a = Nm, — 2.5 x 10"*, Rq = 8.16 
/im, A(j> — 37r/2 and Ax — 0.56 fim. Here, for the ^'^Na 
condensate, the ground state is obtained by the acting of 
C/(0, 7r/2, 7r/4) on the ground state that we obtain from 
direct numerical calculation. 

In this case, the coexistence state structure decays into 
vortex and anti-vortex pairs at an earlier time for the 
^^Rb condensate. As it is shown in Fig. [TOl at about 
t ~ 123.22 ms, four vortex and anti-vortex pairs arrang- 
ing themselves in a "-f " configuration have been formed 
in the tpo component of the condensate while vortex pairs 
are emerging in the other two components. The vortex 
and anti-vortex pairs afterwards arrange themselves in 
-I— like and x-like configurations alternatively. Following 
the vortex configuration in the ipo component, the vor- 
tex and anti-vortex pairs take the configuration of 
in -01 and tp-i components at t = 164.28 ms, while they 
take the configuration of "x" in the tpo component. At 
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FIG. 10: (Color online) The density distributions IV^il^, 
and I'f/'-ip of *^Rb condensate at chosen time points with 
Nm = 2.5 X 10*, Ro = 8.16 urn, Afj] = 37r/2 and Ax = 0.56 
/xm, showing the decay of the ring dark solitons and motion 
of the vortex and anti-vortex pairs, (a)-(c), (d)-(f) and (g)-(i) 
are density configurations at t = 123.22 ms, t = 164.28 ms, 
and t = 189.10 ms, respectively. 



t = 189.10 ms, the vortex and anti-vortex pairs in ipi 
and components catch up with that in the com- 
ponent, and evolve into x-like configurations. However, 
this kind of vortex pairs oscillations do not occur in the 
^^Na condensate with a time span of 378.64 ms, at which 
the ring dark solitons breaks into four pieces in the ipQ 
component as shown in Fig. [II](a)-(c). The insets to the 
top-right of the density distributions show the phase dis- 
tributions, showing the absence of the vortex and anti- 
vortex pairs. The coexistence state composed of ring 
dark solitons coming from different components in ^^Na 
condensate enjoy a longer lifetime than that in the *^Rb 
condensate. We ascribe this phenomena to the effect 
of the spin-dependent interaction parameter C2, which 
is positive for ^'^Na and negative for ^^Rb. 

Taking the same parameters in the *^Rb condensate 
but with a positive as (o^s — 0.038a„ is chosen in our 
numerical calculation), we find that the lifetime of the 
coexistence state composed of ring dark solitons is pro- 
longed remarkably, as show in Fig. [TT] (d)-(f). In the 
real ^^Rb condensate with parameters in Fig. [TUJ the 
coexistence sate breaks into vortex and anti- vortex pairs 
totally at 123.22 ms (Fig. [in](a)-(c)). However, the state 
can last until t = 160.86 ms as shown in Fig. [Til (d)-(f). 
Different spin mixing dynamics which is determined by 
the spin-dependent interaction term are found between 
^■^Na and ^^Rb condensates as shown in Fig. [T^l and 
this should be responsible for the different lifetimes of 
the ring dark solitons between ^"^Na and ^''Rb conden- 
sates. 

So ring dark solitons of the ferromagnetic condensate 
incline to decay into vortex and anti-vortex pairs much 
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FIG. 11: (Color online) (a)-(c) are the density distributions 
IV'il^, liAoP and IV'-il^ of ^^Na condensate at t = 378.64 with 
NNa = 2.5 X lO"*, Ro = 8.16 /xm, A(j) = 3tt/2 and Ax = 0.56 
fj,m. The insets on the top right shows the phase distribution 
of the condensate, (d)-(f) are density distributions at t = 
160.00 ms with the same parameters as in Fig. 1101 except for 
as = 0.038an. The ground sate used have been rotated by 
f/(0, tt/2, it/4:) in the spin space. 



earlier than that of the antiferromagnetic condensate. As 
it is known, modulational instability is possible in ferro- 
magnetic spinor condensates [1^ [3l| . And the instabil- 
ity would cause the breakdown of exotic structures such 
as ring dark solitons. Calculations in the present paper 
show that such a worry is not necessary. The modula- 
tional instability has a close dependence on the phase dis- 
tribution of the state, as found in [U . By simulating the 
experimental generation of ring dark solitons, the state in 
our present setup, which is composed of interdependent 
ring dark solitons, is not a stationary state. So the con- 
dition for the modulational instability in [3l| is not 
well satisfied. We can understand this point by referring 
to the spin-mixing effect. States that would suffer mod- 
ulational instability found in [2^ [sj] do not involve the 
spin-mixing [ssj between components of the condensate. 
In our system, spin-mixing occurs at the very beginning 
as shown in Fig. I12[ which exhibits the spin populations 
versus evolution time (spin-mixing in other parameter 
region has also been found), indicating the deviation of 
the state of ring dark solitons from the modulationally 
instable states obtained in [23, ^ij • 



V. CONCLUSIONS 

In conclusion, the generation and dynamics of the ring 
dark solitons in both spin-1 ^'^Na and ®^Rb condensates 
are explored by simulating the phase engineering tech- 
nique numerically. We find that the characteristics of the 
far-off resonant laser pulse, such as the Stark-shift poten- 
tial width (revealed by Ax) and the duration of the beam 
St (revealed by Ac/)), affect the dynamical properties of 
the condensate dramatically, including the lifetime and 
the decay profiles of the ring dark solitons. We choose 
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FIG. 12: (Color online) Spin populations versus evolution 
time of ^^Na (upper panel) and ^'^Rb (bottom panel) conden- 
sates with NNa = Nm = 2.5 x lO", Ro = 8.16 /xm, Atf) = 3/2tt 
and A.X = 0.56 fim. For the ^^Na condensate, a rotated initial 
ground state is used. 



optimal parameters to generate nearly black ring soli- 
tons. If only one ring dark soliton in one component of 
the condensate is generated, ring dark solitons in other 
components are induced, and they compose a coexistence 
state, which exhibits dynamical oscillations for hundreds 
of milliseconds. The oscillation of the coexistence state 
can live hundreds of milliseconds because of the mutual 



filling and spin mixing of different components. The co- 
existence state is more likely to be stabilized in presence 
of all the three components of the condensate. The co- 
existence state can enjoy a longer lifetime in the ^'^Na 
condensate than that in the ^^Rb condensate. How- 
ever, if the coexistence structure is induced by an ob- 
viously off-centered ring dark soliton, it suffers an earlier 
collapse with the appearance of vortex and anti-vortex 
pairs. This work may help to realize manipulation of 
ring dark solitons in high-dimensional multicomponent 
BECs. Because of the long lifetime of the ring structures 
(more than 800 ms for ®^Rb), we hope that our results 
can stimulate investigation of high-dimensional ring dark 
solitons in future spinor BEC experiments. In fact, the 
density engineering and phase engineering can be utilized 
together to generate more perfect ring dark structures, 
whose lifetime should be much longer. Besides, the FR 
technique is waiting to be used to manipulate ring dark 
solitons in spinor BEC in the future. 
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